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.^ We study conductance fluctuations (CF) and the sensitivity of the conductance to the mo- 

(N 

f~^ tion of a single scatterer in two-dimensional massless Dirac systems. Our extensive numerical 

Q study finds limits to the predicted universal value of CF. We find that CF are suppressed for 

m 

y—{ ballistic systems near the Dirac point and approach the universal value at sufficiently strong 

. |_( disorder. The conductance of massless Dirac fermions is sensitive to the motion of a single 

Vh scatterer. CF of order e /h result from the motion of a single impurity by a distance com- 

parable to the Fermi wavelength. This result applies to graphene systems with a broad range 

of impurity strength and concentration while the dependence on the Fermi wavelength can be 
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explored via gate voltages. Our prediction can be tested by comparing graphene samples with 
varying amounts of disorder and can be used to understand interference effects in mesoscopic 
graphene devices. 

keywords: graphene • transport property • conductance fluctuations • electronic transport in nanoscale 
materials and structures ■ nanoelectronic devices 

The two-dimensional (2D) Dirac equation is relevant in elucidating the electronic and transport 
properties of recently discovered materials such as graphene and topological insulators. Graphenej'^S 
a one- atom- thick allotrope of carbon, and the topological surface states of materials with a bulk 
gap such as Bii__^Sbx crystals'^) share unusual Dirac-like electronic structure providing an en- 
thralling test bed for new physics and inevitable future applications based on quantum interference 
effects.S'^ Given the Dirac nature of the electronic spectrum, the quasiparticles in graphene prop- 
agate as massless relativistic particles, making graphene a qualitatively different material from 
conventional conductors. At low temperatures, the quantum transport of electrons becomes co- 
herent and leads to quantum interference corrections to the conductance. The amplitude of these 
corrections follows a universal scaling with value dG ~ e^/h receiving the name of universal con- 
ductance fluctuations (UCF).^^ The conductance fluctuations (CF) take place when a coherent 
electron wave scatters repeatedly while traversing a disordered conductor. The wave follows all 
possible paths through the sample and different paths interfere with each other, giving rise to CF 
that are independent of the sample size and the degree of disorder. Due to the CF, the conductivity 



is sensitive to changes in the configuration of the impurity scatterers. Theoretical studiesl^CSIhave 
predicted that the full UCF effect is obtained by moving a single scatterer a distance comparable 
to half the Fermi wavelength Xf. In gate-doped graphene the gate voltage controls the average 
charge density of the device and given the relationship Xf = 2(;r/|n|)'''^, the gate also presents an 
exquisite control over the Fermi wavelength. The Fermi wavelength diverges as the carrier density 
approaches the neutrality (Dirac) point, making the low-density limit particularly interesting given 
that experimentally the conductance does not go to zero and disagreements between theory and 
experiments regarding the minimum value of the conductivity exist. 



Recent experimental and theoretical studies of the CF in graphene have been reported.S^'^^n^'^ 
Berezovsky and coUeaguesQSEQ] demonstrated unambiguously the sensitivity of the CF to the mo- 
tion of a single scatterer. Placing a charged scanning probe microscope tvfi^ near the graphene 
sample creates an image charge that acts as a movable scatterer. By mapping the CF versus scat- 
terer position, it was found that the UCF decorrelate when the induced scatterer is displaced by a 
distance comparable to half the Fermi wavelength.QSl 

In this article we calculate the zero temperature CF for quasiparticles obeying the massless 2D 
Dirac equation in the presence of disorder, assuming elastic scattering by fixed scatterers and with 
no sources of inelastic scattering. The motivation for the study of this system is threefold: First, 
the massless Dirac Hamiltonian model is of relevance to transport in graphene and the surface 
states of topological insulators. Second, we probe the sensitivity of the conductance to changes 
in the impurity strengths over a range of carrier densities and system sizes that are accessible by 
current experiments. Finally, we calculate the CF due to the motion of a single scatter. This last 
result demonstrates the sensitivity of CF to the motion of an ionized charged impurity and shows a 
dependence on the carrier density that applies to all but the most disordered graphene devices. 

Results and Discussion 

Model. We study the following 2D Hamiltonian: 

H = -inv{a^d^ + aydy) + V{r)ao (1) 

with V being the velocity of the Dirac fermions, (T, the 2x2 Pauli matrices, and V (r) a (pseudo)spin- 
independent potential. The spinor wave function applies to graphene near energies close to the 
Dirac (neutrality) points with each component pertaining to one of the two atoms in the unit cell. 
Our study is based on a transfer matrix approach to calculate the conductance of Dirac fermions,'22l 
described in the methods section. The conductance is modeled by calculating the reflection r and 
transmission t matrices in realistic devices. Although the use of the Dirac equation does not take 



into account the anisotropy of the Dirac cones (trigonal warping) in graphene, the resulting trigonal 
distortion occurs at doping values considerably larger than those presented here. 

Experiments have observed an electrostatic landscape with naturally occurring variations in the 
carrier density .'^^ffiSI jj^ (^j^^^ study we consider such a potential arising from the random charged 
impurities located either in the substrate or above the graphene plane.'^^ES This type of disorder 
profile is both smooth and slowly varying on the atomic scale, suppressing intervalley scattering. 
The disorder model used in our calculations neglects scattering from short-range atomic defects 
and ripples. This approach is justified by the long-range (compared to the atomic lattice) nature of 
the Coulomb-like charged impurities. 

After choosing the average charge density hq we introduce local variations in the charge den- 
sity either through randomness of the on-site carrier density at each lattice point rjj or from ion- 
ized impurities near the graphene plane. In the first case, the contribution of the spatially vary- 
ing on-site carrier density is distributed uniformly with width An [^(r) = n{rj) = no + rjj with 
—An < rjj < An]. In the latter procedure, a fraction n, of the lattice sites are randomly chosen 
and an ionized impurity center is located a distance away. We solve for the charge density in- 
duced from this set of ideally screened Coulombic scatterers. Using the relationship between 
the Fermi energy and carrier density we obtain the local scattering potential that is used in [T] 



V{rj) = fiv ^%n{n{r j)) ^J Ttn{r j) . Note that the square root means that the contributions to the po- 
tential do not add arithmetically. In this simple model the random on-site component dominates 
at low carrier density and its contribution diminishes as the overall charge density increases. Al- 
though our approach overestimates screening at the interface of hole-rich/electron-rich regions, the 
charge of the scatterers is chosen to yield an rms charge density in agreement with the observed 
charged puddles. l^^^ A more rigorous (quantum) treatment!^ of the density response is limited to 
system sizes significantly smaller than those we have investigated in the present work. 
Conductance. The conductance is obtained from the Landauer formula. 
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Figure 1: Conductance fluctuations of Dirac fennions. The CF are plotted in terms of the standard 
deviation of G (5G in e^/h) versus carrier density n (in cm^^) for several values of the fluctuations 
in the carrier density width An (from bottom to top An = 1.0,2.0,3.0, and 4.0 x lO^^cm"^) and 
system sizes (a) L = 170, (b) L = 350, and (c) L = 670 nm. (d) Transmission per mode at the Dirac 
point for different values of An. After the first mode, all other modes come in degenerate pairs 
with only one shown in the plot, (e) Scaling of the CF with system size for different values of An. 
The length of the sample is varied with a fixed W/L ratio of 3. 



where Go is the quantum of conductance {4e^/h for graphene due to spin and valley degeneracy), 
N is the total number of transverse modes in the sample, and r„ are the transmission eigenvalues 
obtained from the diagonalization of the matrix product tt\ The transmission matrix t is obtained 
using the transfer matrix approach. The details of the calculation are described in the methods 
section. 

Conductance fluctuations. [T](a-c) presents the CF, with dG defined as the standard deviation 
of the conductance, as a function of carrier density in the no decoherence limit. Measurements 
of the magnitude of the CF at the Dirac point have yielded mixed results. In weakly disordered 
bilayer and trilayer graphene devices the UCF are suppressed near the Dirac point.^^ Although 
the presence of more than one layer of graphene substantially changes the electronic properties 
of the device,'^ possible explanations for the suppression of the fluctuations'^ involved a different 



mechanism for quantum interference of edge states, which dominate the conductance of the low- 
density samples. A reduction of the amplitude of the fluctuations was seen in monolayer graphene 
devices'^2112 contrasting with a study^ in which the amplitude of the fluctuations is larger at the 
charge neutrality point for both monolayer and bilayer graphene devices. Our calculations, cov- 
ering a broad range of impurity strengths and concentrations, finds that the amplitude of the con- 
ductance fluctuations at the Dirac point for sufficiently strong fluctuations in the carrier density 
An > 2.0 X lO^^cnT^ exhibits a peak. In contrast, the CF decrease at exactly the Dirac point for 
weak on-site density fluctuations. 

The suppression of the CF in weakly disordered systems is due to the reduced size of the 
fluctuations in the cleaner systems coupled to a reduced number of modes open to conduct's^ 
[T] (d) shows the transmission for the first modes at the neutrality point for several values of the 
on-site carrier density averaged over several disorder configurations in L = 350 nm systems. All 
systems, independent of the amount of disorder present, have at one mode completely open. In 
the clean systems, the second mode is only partially open {T = 0.06). Disorder opens up the 
subsequent modes for transport as can be seen in [T] (d) where error bars show the fluctuation for 
each of the modes. Assuming that the CF are proportional to the fluctuations in each of the modes, 
the suppression of the CF is due to the smaller fluctuations in the cleaner systems coupled to the 
reduced number of modes open to transport. 

In [T] (e) we plot the averaged CF as a function of system size. For systems with a fixed W/L 
ratio of 3, we vary the length from 170 nm to 670 nm. As the strength of the on-site carrier density 
increases, the value of 5G converges to 2. 17 ± 0. 15 e^ /h, close to the value found analytically for 
doped graphene of 2.36 e^ /h in the W/L 3> 1 aspect ratio limit.^^^^^ The value is higher than the 
UCF value {e^ /h) due to the absence of intervalley scattering and trigonal warping, both absent in 
our model. 

Sensitivity to motion of a single impurity. Similarly, we compare the CF from different im- 
purity configurations to the CF induced by the motion of a single impurity. The method developed 
by Tworzydlo et al.^ is suitable to revisit a landmark study carried out by Feng et al. for metallic 
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Figure 2: Conductance fluctuations due to the interchange of two on-site densities. The rms devi- 
ation of the conductance upon interchange of the local charge density of two sites (rfj) for several 
values of the fluctuations in the carrier density with width An (from bottom to top An= 1.0,2.0,3.0, 
and 4.0 x lO^^cm^^) and system sizes (a) L = 170, (b) L = 350, and (c) L = 670 nm. 



systems.f^We obtain the conductance of massless Dirac fermions for a discrete lattice model with 
random on-site charge densities (similar to the Anderson model) and calculate the rms deviation 
of this conductance dGi after the interchange of on-site charge densities of just one pair of sites. 
In [2] we plot dGi as a function of carrier density for several system sizes averaged over different 
realizations of the random densities and for several widths An over which these on-site values are 
chosen. While an increase in the strength of the on-site carrier density produces an increase in the 
rms value of CF, we found 5Gi is independent of sample size. Our calculation chooses two sites 
to interchange at random and as such their separation can be small or their on-site values similar. 
Thus, although the average value of 5Gi is lower than the average value of dG, in general dGi is 
bounded by dG. 

For weak disorder and small system size (An = 1 .0 x 10^^ cm^^ and L = 170 nm), CF for both 
the complete configuration change (lower curve in[T](a)) and upon interchange of two sites (lower 
curve in |2] (a)) exhibit oscillations as a function of the carrier density. In such systems ballistic 
transport dominates, i.e., the transport mean free path is larger than the system size / > L. The 
oscillations in the CF as a function of density are caused by multiple reflections at the ends of the 
sample (Fabry-Perot resonances) where the enhancements are due to multiple visits .1^^^ 

Finally, we consider the change in conductance induced by the motion of a single charged 



impurity a distance 5r. A fraction of lattice sites are randomly selected and charged impurities 
are placed above their positions (ri,?"2,---5'"A?,) inducing a charge density landscape such as the 
one presented in [3] (a). Such an experiment was undertaken by Berezovsky et al. where the 
device conductance was measured as a function of SPM tip position, and the movable scatterer 
is created by the SPM tip.^^l We present such conductance maps obtained from raster scanning 
one of the impurities in|3](b). Maps are presented for two different carrier densities and show 
that the lateral size of the fluctuations depends on the carrier density. It was shown that the length 
Ir that a strong scattering center needs to be shifted to decorrelate the CF (i.e., to change the 
conductance by ~ e^/h) is approximately half the Fermi wavelength.-^ This length was obtained 
from the autocorrelation function of the conductance maps and studied as a function of the carrier 
density.™ In this new study, we calculate the conductance change for the motion of one impurity. 
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over a range of system sizes and carrier densities for several impurity densities, accessible to 
current experiments. 

For diffusive metallic systems 5Gi is given by*^ 
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where Q. is the volume and d the dimensionality. In 2D systems, the L/l factor drops out and 
Q./Nfi corresponds to the strength of the impurities. |4] predicts that if kpdr > 1 the change in 
conductance from the motion of a single impurity will be comparable to the complete change of 
the impurity configuration in a sample. In[3](c) we show the rms fluctuations of the conductance 
as a function of carrier density and distance moved by the impurity. The dashed line in the plot 
indicates the length l^ = O.A6Xf separating the regions at which motion of a single impurity will 
impact the CF as if the impurity configuration was completely different. As evident in[3](c), these 
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Figure 3: Evidence of the conductance sensitivity of Dirac fermions. (a) Typical carrier density 
fluctuations resulting from random charged impurities localized in the substrate or above the plane 
in a graphene device. The concentration of charged impurities is 20 % of sites, and are located 
1 to 10 nm away from the graphene plane with a charge of ±2 e, yielding density fluctuations 
n{r) ~ 1.6 X lO^^cm^^ for system size W = 1.005 iim and L = 670 nm. (b) Conductance maps, the 
device conductance versus the position of the movable scatterer, at two different carrier densities 
simulated by raster scanning a single charged impurity over a 400 x 400 nnP' area of the sample. 
The maps display spatial CF dG\ ~ e^ /h with the lateral size of the fluctuations features depending 
on the carrier density, (c) Plot of the conductance change 5G\ versus the charge density n and 
the distance dr that a single ionized impurity moved, averaged over 2000 impurity ensembles. 
The black dashed line is presented as a guide indicating the length Ir where UCF theory predicts 
that the motion of a single impurity will induce a change in the CF comparable to that of a new 
configuration of impurities. 



computed results indicate a deviation from conventional UCF theoryO^lfor 2D electron gas of Dirac 
fermions. In connection with previous studies, ^21 the Ir length obtained from the autocorrelation 
function will saturate at certain values of the carrier density. One of the main findings of this letter 



is that the trend presented in[3](c) applies to samples where the carrier-density fluctuations ranges 
from 10^ to 3 X lO^^cm ^, i.e., for all but the dirtiest graphene devices. 

Conclusion 

In conclusion, we have focused on the CF of Dirac fermions and the impact of the motion of a sin- 
gle impurity. Our results predict that the CF in 2D Dirac systems is dependent on the strength of the 
disorder near the neutrality point in ballistic systems but independent of the strength of the disorder 
for doped graphene. For strong enough disorder, as is the case of graphene on Si02 substrates, the 
CF will not depend strongly on carrier density but are enhanced at the Dirac point. Studying CF 
on suspended graphene'221 and graphene deposited on hexagonal boron nitride substrateslSScan test 
these predictions. Consistent with theoretical predictions, the change in conductance caused by the 
motion of a single impurity (covering a small area of the sample) is significant when the distance 
moved is of the order of the Fermi wavelength. 

Methods 

Transfer matrix approach. We sketch the transfer matrix procedure used, for further details we 
refer to Tworzydlo et al. ^^ The system is discretized into a lattice and the difference equations 
are solved without violating symplectic symmetry and current conservation while avoiding the 
fermion doubling problem.221 We calculate the conductance in a strip geometry, discretizing the 
sample using a square lattice. The longitudinal direction extends from x = io x = L and the 
transverse direction from y = Q\.oy = W , where L and W are the length and width of the sample. 
Periodic boundary conditions are used in the transverse direction. The transfer matrix reads 
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where ^',„ is a vector containing the values for the wavefunction ^'{x,y) at jc = mao, with m being 
an integer and gq = 10 nm the lattice spacing. Semi-infinite metallic leads are attached to the 
strip at its ends (x = and x = L). The metal contacts are ballistic leads in which all modes are 
conducting.221 Each incoming mode on a lead is propagated to the other lead using the transfer 
matrix. The N transverse modes in the sample are either propagating modes (^/ or evanescent 
modes Xi (modes that decay for large positive or negative values of x). An incoming wavefunction 
in mode Iq starting on the left side of the sample is composed of incoming, reflected, and evanescent 
modes 

^/o (^ = 0) = (^+0 + £ r;,/„ ^- + £ a;,z„zr (6) 

at the X = L edge of the sample the wavefunction is the sum of the transmitted and evanescent 
waves given by 

'^/o(^ = ^) =E^/,/o<^/+ + E</„Zr (7) 

where the label + corresponds to right moving and — to left moving modes. The reflection r/ /^ 
and transmission coefficient ?/ /^ are obtained from the transfer matrix relation 

%^{x = L)=^^i^{x = 0) (8) 

and elimination of the a and a' coefficients. Once this is done for all possible modes the reflection 
r and transmission t matrices are composed. Similarly, repeating the procedure for modes propa- 
gating from the right to the left edge of the sample results in the r' and t' matrices.!^ The conduc- 
tance is calculated using the Landauer formula, [2[ by summing over the transmission eigenvalues 
T„ obtained from the diagonalization of the matrix product tt^ . 
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